Low-decoherence flux qubit 



J. Q. You, 1,2 Xuedong Hu, 3 S. Ashhab, 2 and Franco Nori 2,4 

Department of Physics and Surface Physics Laboratory (National Key Laboratory), Fudan University, Shanghai 200433, China 
Frontier Research System, The Institute of Physical and Chemical Research (RIKEN), Wako-shi 351-0198, Japan 
3 Department of Physics, University at Buffalo, SUNY, Buffalo, NY 14260-1500, USA 
■* Center for Theoretical Physics, Physics Department, Center for the Study of Complex Systems, 
University of Michigan, Ann Arbor, MI 48109-1040, USA 
(Dated: February 6, 2008) 

A flux qubit can have a relatively long decoherence time at the degeneracy point, but away 
from this point the decoherence time is greatly reduced by dephasing. This limits the practical 
applications of flux qubits. Here we propose a new qubit design modified from the commonly 
used flux qubit by introducing an additional capacitor shunted in parallel to the smaller Josephson 
junction (JJ) in the loop. Our results show that the effects of noise can be considerably suppressed, 
particularly away from the degeneracy point, by both reducing the coupling energy of the J J and 
increasing the shunt capacitance. This shunt capacitance provides a novel way to improve the qubit. 



Superconducting quantum circuits based on Joseph- 
son junctions (JJs) are promising candidates of qubits for 
scalable quantum computing (see, e.g., [H). Like other 
types of superconducting qubits, flux qubits have been 
shown to have quantum coherent properties (see, e.g., 
@, S, L3, 0, S 0, 0)- A recent experiment 0] showed that 
this qubit has a long decoherence time T 2 (~ 120 ns) 
at the degeneracy point; this T 2 can become as long 
as ~ 4 /us by means of spin-echo techniques. However, 
even slightly away from the degeneracy point, the de- 
coherence time is drastically reduced. This sensitivity 
to flux bias considerably limits the applications both for 
flux qubits for quantum computing, and also when per- 
forming quantum-optics and atomic-physics experiments 
on microelectronic chips with the qubit as an artificial 
atom. 

Typically, JJ circuits have two energy scales: the 
charging energy E c of the J J, and the Josephson cou- 
pling energy Ej of the junction. Ordinarily, a flux qubit 
works in the phase regime with Ej/E c 3> 1, where its 
decoherence is dominated by flux fluctuations. For the 
widely used three-junction flux qubit design 0, H, 0, H, Hi , 
in addition to two identical JJs with coupling energy Ej 
and charging energy E c , a third J J, which has an area 
smaller by a factor a ~ 0.7, is employed to properly ad- 
just the qubit spectrum. Charge fluctuations can affect 
the decoherence of this flux qubit via the smaller junc- 
tion. 

Here we search for an improved design for flux qubits. 
We show that reducing the ratio Ej/E c suppresses the 
effects of flux noise, although charge noise becomes in- 
creasingly important. Reducing a further suppresses the 
effects of flux noise and considerably improves the deco- 
herence properties away from the degeneracy point. As 
the effect of flux noise has been largely suppressed, charge 
noise would now be the dominant source of decoherence. 
It mainly comes from the charge fluctuations on the two 
islands separated by the smaller JJ and affects the qubit 
mainly through relaxation. We thus propose an improved 



flux qubit by introducing a large capacitor that shunts in 
parallel to the smaller J J. This shunt capacitance sup- 
presses the effects of the dominant charge noise in the 
two islands separated by the smaller JJ by reducing the 
charging energy. Our results reveal that using a larger 
shunt capacitor allows reducing both Ej/E c and a to 
considerably suppress the effects of both flux and charge 
noises, particularly away from the degeneracy point. In 
essence, our method reduces the couplings of the flux 
qubit to the two types of noise. It provides a promising 
approach for lowering the decoherence of J J qubits. 

We consider, as a typical example, a modified version 
of the usual three-junction flux qubit [1, H, 0, H, [H] ; it con- 
sists of a superconducting loop containing three JJs and 
pierced by a magnetic flux. The new ingredient, which 
drastically improves the qubit, is an additional capaci- 
tance C s shunted in parallel to the smaller J J (Fig. [T]). 
Note that our approach also applies to flux qubits with 
any number of junctions, e.g., the four-junction design, 
by shunting a capacitor to the smaller junction. As 
shown in Fig. [TJ the three JJs divide the superconduct- 
ing loop into three islands, denoted by a, b, and c. When 
the environment-induced charges on the three islands are 
taken into account, the Hamiltonian of the system is H = 
E p (n p - SN a ) 2 + E m [n m - (SN h - SN C )] 2 + U(cp p ,<p m ), 
where — —id/dtpk {k — P,m)i E p — 2E C , and E m = 
E p /(1 + 2f3), with E c = e 2 /2C, 7 , and (3 = a + C s /Cj. 
Here E p is the charging energy of the island a, E m is 
the effective charging energy related to islands b and c, 
and <EVj (i — a, b and c) are the environment-induced 
charges (in units of 2e) on the islands. The poten- 
tial energy is U((p p ,ip m ) — 2Ej(l — cos ip p cos <p m ) + 
aEj[l - cos(27r/ + 2ip m )], where ip p = (ipt + (p 2 )/2, 
and tp m — (ipx — ip 2 )/2. The reduced flux / is given 
by / = fe + Sf, with / e = $ e /$ and Sf = <5$/$ , where 
$ e is the externally applied magnetic flux in the loop, 
the flux fluctuations, and $ = h/2e. 

We study the environment-induced effects on the qubit 
by expanding the Hamiltonian H into a series and trun- 
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FIG. 1: (Color online) (a) Schematic diagram of a low- 
decoherence flux qubit. The loop is pierced by an externally 
applied magnetic flux <J> e and interrupted by three Joseph- 
son junctions (JJs) that divide the loop into three islands 
denoted by a (blue), b (red) and c (red). The two identical 
JJs have coupling energy Ej and capacitance Cj, while the 
third (smaller) JJ has coupling energy aEj and capacitance 
aCj, where 0.5 < a < 1. Here a capacitance C s is shunted 
in parallel to the smaller J J, so as to reduce the charging en- 
ergy related to islands b and c. The shunt capacitor can be 
fabricated by depositing a metallic plate below or above the 
smaller junction, with a dielectric layer between the metallic 
plate and b as well as between the plate and c. 
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eating it to keep the interaction terms to lowest order. 
Including the environmental Hamiltonian H cnv , we can 
write the total Hamiltonian H t as H t = H q + H cnv + Hi nt , 
where H q = H(8Ni = 0,Sf = 0) is the Hamiltonian of 
the qubit and the interaction Hamiltonian takes the form 
Hut = -2E p n p SN a -2E m n m {8N b -5N c )-I $ Sf, where 
I = -a I c sin(27r/ e + 2cp m ) with I c = 2t:Ej/<S> . 

To study the coherence properties of this flux qubit, we 
project the total Hamiltonian onto the subspace spanned 
by the qubit eigenstates |0) and |1) with eigenenergies 
£0 and Ei, the two lowest energy levels of the quantum 
device. Now the flux qubit Hamiltonian is reduced to 
H q = \ea z , with e = E x - E and a z = - |0)(0|; 

while the interaction Hamiltonian is reduced to 

H int = -J2[^X i -(a + Y i + R.c.)], (1) 

i 

with <t+ = |1)(0|. Here the longitudinal couplings Xi, 
i = f , a, b and c, are Xi(t) = AiSXi(t), with Af = 
±$o«l|/|l}-<0|/|0)), A a = E p ((l\n p \l)-(0\n p \0)), 
and A b = -A c = E m «l|ra m |l} - (0|n ro |0)). The 
transverse couplings Yi are Yi(t) — Bi5Xi(t), where 
B f = $ (1|/|0>, B a = 2E p (l\n p \0), and B b = -B c = 
2E m (l\n m \0). The fluctuations are SXf = Sf for the flux 
noise and SXi = 8Ni (i = a, b, and c) for the charge noises 
related to the three islands. The longitudinal coupling 
term a z Xi leads to pure dephasing between the qubit 
states, while the transverse coupling term, a+Yi + H.c, 
leads to relaxation. One way to suppress decoherence 
from both pure dephasing and relaxation is to reduce 
the longitudinal and transverse couplings by decreasing 
\Ai\ and \Bi\. This general method of decoherence sup- 
pression applies irrespective of the particular behavior of 



FIG. 2: (Color online) (a)-(e): Energy levels of the flux qubit 
versus the reduced magnetic flux f e = <E> e /<I>o for a = f3 = 0.8 
and Ej/E c = (a) 60, (b) 35 and (c) 20 as well as a = (3 = 0.6 
and Ej/E c — (d) 60 and (e) 35. Here only the lowest four 
energy levels are shown, (f)-(j): Coupling strength (in units 
of Ej) versus the reduced flux f e , which correspond to (a)- 
(e), respectively. In this figure, and the following one, other 
coupling strength parameters for charge noise are not shown 
because they are orders of magnitude smaller. 



Aj(t), i.e., whether it is Gaussian or non-Gaussian noise. 

We first study the conventional three-junction flux 
qubit, without the shunt capacitance C s . We take 
a = /3 = 0.8, as in the experiment in Ref. Q. Fig- 
ures [UJa)-[2Ic) show the flux dependence of the energy 
levels. The lowest two levels around the degeneracy point 
/ e = & e /&o = 0.5 are employed as the qubit states. To 
characterize the effects of the flux noise on this qubit, we 
introduce longitudinal and tranverse coupling strengths 
defined by F z = \Af\ 2 and F x = \Bf\ 2 , respectively. 
To compare the contributions of charge noise with flux 
noise, we define the longitudinal and transverse coupling 
strengths C Z i = Ki \Ai\ 2 and C X i = re, |-B;| 2 for the charge 
noise, where i = a, b and c. The coefficient k% character- 
izes the relative contribution of each charge noise and is 
defined as Ki — Si(co)/Sf(u>), i.e., the ratio between the 
power spectra of each charge noise and the flux noise. 
This definition is reasonable as the qubit relaxation rate 
is proportional to both \Bi\ 2 and the power spectrum of 
the noise, as shown below [see Eq. (j2])]. We estimate 
Ki by considering the 1// noise with power spectrum 
Ki/\oj\. Here Ki is determined from experiments. Typi- 
cally, -Kcharge = (0.3 x 10 -3 ) 2 for the charge noise [9( and 
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-Kfiux = 3 x 10~ 12 for the flux noise 0|. For simplicity, the 
same -Kcharge is used for the three charge noises related 
to the islands a, b and c, so that Ki = -Kcharge / Kfiux = K - 

In Figs. Oi)-lHh), we show the flux dependence of the 
coupling strengths F z , F x , and C x . At f e = 0.5, F z falls 
to zero, while F x rises to its peak. This implies that, at 
the degeneracy point f e = 0.5, the first-order pure de- 
phasing due to flux noise disappears, and the flux qubit 
decoherence is dominated by relaxation. However, for a 
large Ej/E c [see, e.g., Fig. HJf)], pure dephasing domi- 
nates when f e is slightly away from the degeneracy point. 
For decreasing Ej/E c , the valley of F z around f e = 0.5 
becomes broader, while the peak of F x becomes less sharp 
and its height is gradually reduced. In Figs. E^f )-Hpi) we 
also show the dominant C x b = C xc = C x curves due 
to charge fluctuations on the smaller islands b and c. 
The quantity C x characterizes qubit relaxation induced 
by charge noise and should be compared to the F x curves 
in the figures. Notice that when Ej/E c decreases, charge 
noise plays an increasingly important role and eventually 
is more important than flux noise in terms of relaxation 
[compare the dashed and dotted curves in Fig.^h)]. We 
note here that as shown in 0, the energy levels of the 
flux qubit are very flat versus the offset charges. This 
insensitivity of the energy splittings to the charge fluc- 
tuations implies a very weak pure dephasing caused by 
charge noise. Indeed, we numerically calculated the lon- 
gitudinal coupling strengths C za and C z b — C zc = C z for 
our system, and found that these quantities are orders of 
magnitude smaller than C x . These results further reveal 
that the charge-noise-induced dephasing is weak. 

Now we reduce the size of the smaller J J in the conven- 
tional 3-JJ flux qubit to a = = 0.6. We also show the 
flux dependence of the energy levels [Figs.[2]Jd) andEte)] 
and the coupling strengths [Figs.^i) andGDj)]. Here the 
coupling strengths F z and F x are reduced and become 
flatter around / e = 0.5. This reduction of the flux-noise 
effects is due to the decrease of the circulating current 
/ with a. However, the coupling strength C x becomes 
much larger than F z and F x because the effects of charge 
noise are strengthened due to the increase of both charg- 
ing energy E m and transition matrix element |(l|n m |0)|. 
The results here reveal that the decoherence of the flux 
qubit is sensitive to the values of a and (3 of the small 
J J (here a = (3). Moreover, the coupling strength C x 
increases rapidly with decreasing Ej /E c [see the dashed 
curves in Figs. [S^i) and[2]j)]. In short, decreasing a re- 
duces the coupling strength between the flux qubit and 
the flux environment, and makes the coupling strength 
less sensitive to the flux bias so that pure dephasing be- 
comes less important near the degeneracy point. How- 
ever, decreasing a also leads to a dramatic increase in the 
coupling strength between the flux qubit and its charge 
environment, to the degree that it may become the dom- 
inant decoherence channel. 

To achieve an improved flux qubit in which the effects 
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FIG. 3: (Color online) (a) and (b): Energy levels of the flux 
qubit versus the reduced magnetic flux f e . Here, (a) a = 0.8, 
f3 = 1.5 and Ej/E c = 20; (b) a = 0.6, j3 = 4 and Ej/E c = 35. 
(c) and (d): Coupling strength (in units of Ej) versus the 
reduced flux f e , which correspond to (a) and (b), respectively. 

of both charge and flux noises are reduced significantly, 
we shunt a large capacitance in parallel to the smaller 
JJ (see Fig. [T]) so as to decrease the charging energy E m 
while keeping the ratio a small. In Fig. [3] we present 
two examples in which the effects of the charge noise are 
reduced. In both small and flat C x is achieved. 

Also, we show that the coupling strengths F z and F x are 
smaller and flatter in Fig. 02d) than in Fig. |3](c). These 
results indicate that for a suitably chosen Ej/E c ratio, by 
optimally decreasing a and increasing [3 one can reduce 
the coupling of the qubit to both flux and charge noises, 
so that pure dephasing can be considerably reduced in a 
wide region around the degeneracy point f e — 0.5 and the 
relaxation is significantly suppressed. This corresponds 
to an improved flux qubit with low decoherence. 

Note that the parameter a has a lower bound of 0.5 
for the flux qubit; when a < 0.5, the double-well poten- 
tial reverts back to a single- well potential and the circuit 
behaves like a phase qubit. Also, the shunt capacitance 
should have an upper bound when other factors are taken 
into account. For instance, a very large shunt capacitance 
needs a thicker dielectric insulator for fabricating the ex- 
ternal capacitor. In this case, the decoherence originating 
from phonon radiation [Io| and defects in the thicker in- 
sulator [11] may play more important roles. Furthermore, 
when the shunt capacitance increases, the energy gap A 
of the qubit at the degeneracy point narrows down, rais- 
ing the single-qubit operation time h/A. This decrease in 
A is another factor one needs to consider in determining 
the upper bound of the shunt capacitance. One should 
keep in mind here that reducing a increases A. As a re- 
sult, with properly chosen values of a and (3, the gap A 
is not necessarily decreased in the optimized design, as 
we shall show with an example below. 

Finally, we discuss the relation between the decoher- 
ence rate and our defined coupling strengths. When each 
transverse coupling term in Eq. ([T]) is treated as a per- 
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FIG. 4: (Color online) Relaxation and decoherence times (Ti 
and Tb) for three flux qubits with the same noise sources. The 
parameters are Ej/E c — 35, a = ft = (a) 0.8, (b) 0.6, and 
(c) a = 0.6, /3 = 4. Here, as a simple estimation, we ignore the 
frequency dependence of the power spectra around f e = 0.5 
and determine their high- and low-frequency components via 
Eqs. ([2]) and ([3]) by assuming Ti = 0.5 fjs and T2 = 0.8 fis at 
fe = 0.5004 for the qubit in (a). 



turbation, according to the Fermi golden rule, one can 
obtain the (\Bi | dependent) relaxation rate for each noise: 



r 



\Bi\ 2 Si(w 10 ), 



(2) 



where wxo — {E\ — Eq)/%, and the power spectrum is 
denned by S,(uj) = /+~ dt(Xi(t + r)A l (t)) e - iwr . The 

sum of all Ij gives the total relaxation rate Ti = 1/T%, 
where Ti is the relaxation time. 

The longitudinal qubit-environment coupling intro- 
duces a random phase between the qubit eigenstates. 
At time r this random phase is = (1/ti) Jq dtX^t). 
For a Gaussian noise, the dephasing factor 77, (r) is given 
by 0: e-^M = (e lA ^) = exp[-i((A</> i ) 2 }], where the 
brackets (...) denote the quantum statistical average over 
the environment. Because Xi(t) = AiS\i(t), one can 
write the dephasing factor as 



1 



dio Si(ui) 



sin 2 (wT/2) 
2vr(w/2) 2 ' 



(3) 



where tu c is a low-frequency cutoff determined by the 
measurement time. The pure dephasing rate = 1/T V 



is defined by J^iViiT^ 



1. The dephasing factor 



?7,;(r) oc \Ai\ 2 , consistent with reducing pure dephasing 
by decreasing the longitudinal coupling. Moreover, relax- 
ation can also cause damping of the off-diagonal density 
matrix elements. Following the Bloch-Redfield theory 
(see, e.g., [Hj]), the total damping of the off-diagonal ele- 
ments is characterized by a decoherence rate r 2 = l/?2, 

with r 2 = |ri + r v . 

In Fig.[4l we show T\ and T 2 for three hypothetical flux 
qubits calculated using the same noise sources. The qubit 
in Fig. [UJa), with a = (3 = 0.8, is a conventional qubit 
with no shunt capacitance The qubits in Figs. H£b) 
and QJc) have a reduced value of a (0.6 for both fig- 
ures), and a shunt capacitor (with (3 = 4) is added for 



Fig. HJc). As compared with the conventional qubit in 
Fig.[4ja), the decoherence time T 2 in Fig.[4jc) IS larger by 
a factor of 2 at the degeneracy point, and the T 2 peak is 
broader. Furthermore, the qubit in Fig. Etc) has a larger 
gap than the qubit in Fig. Ufa), A w 0.02i5j compared 
to A » 0.013i?j. Combining the increase in both T 2 and 
A, we find that the quality of quantum coherence is im- 
proved by a factor of 3 at the degeneracy point. Away 
from the degeneracy point (e.g., at f e — 0.4988), T 2 can 
be improved by a factor of 3 and the quality of coherence 
is improved by a factor 5 (nearly one order of magnitude). 
Also, at f e ~ 0.5, T 2 is reduced by about one order of 
magnitude if the shunt capacitor is removed [compar- 
ing Figs. BJc) and[UJb)]. These results further show the 
important role of the shunt capacitance in achieving a 
low-decoherence flux qubit. 

In conclusion, we have proposed a new qubit design 
modified from the commonly used flux qubit. The qubit 
decoherence is reduced by shunting the small JJ with an 
additional capacitor. We show that by increasing the 
shunt capacitance and reducing the coupling energy of 
the J J, the effects of both charge and flux noises are con- 
siderably suppressed. Recently, a shunt capacitor was 
used to improve the performance of phase qubits [HI]. 
In that case the motivation for adding the shunt capac- 
itor was quite different from ours; they used a smaller 
junction so as to reduce the number of two-level systems 
in the junction but decreased the charging energy of the 
junction with a shunt capacitor in order to push the qubit 
back into the phase regime. However, in our case the ef- 
fects of the noise are suppressed even though we assume 
that the noise source remains unchanged. 
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